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ABSTRACT 

Using radiative magneto-hydro dynamic simulations of the magnetised solar photosphere and de¬ 
tailed spectro-polarimetric diagnostics with the Fel 6301.5A and 6302.5A photospheric lines in the 
local thermodynamic equilibrium approximation, we model active solar granulation as if it was ob¬ 
served at the solar limb. We analyse general properties of the radiation across the solar limb, such 
as the continuum and the line core limb darkening and the granulation contrast. We demonstrate 
the presence of profiles with both emission and absorption features at the simulated solar limb, and 
pure emission profiles above the limb. These profiles are associated with the regions of strong linear 
polarisation of the emergent radiation, indicating the influence of the intergranular magnetic fields 
on the line formation. We analyse physical origins of the emission wings in the Stokes profiles at the 
limb, and demonstrate that these features are produced by localised heating and torsional motions in 
the intergranular magnetic flux concentrations. 

Subject headings: Sun: Photosphere — Sun: Surface magnetism — Plasmas — Magnetohydrodynam¬ 
ics (MHD) 


1. INTRODUCTION 

The solar photosphere at the limb remains one of the 
largely unexplored solar regions due to the lack of bright 
and extended emissive features, such as spicules in the 
chromosphere, and sharp limb darkening, associated with 
low light conditions and high instrumental noise levels 
in the observations of this region. On the other hand, 
the solar limb has potential for the diagnostics of paral¬ 
lel to the solar surface photospheric flows and magnetic 
fields. Together with solar disk observations, it can pro¬ 
vide new and more detailed information on the geometry 
of the flow and magnetic field structures in the solar pho¬ 
tosphere, on the details of transition between the photo¬ 
sphere and the chromosphere, as well as on the origins 
of chromospheric features. 

Groundbreaking spectro-pola rimetric observation of 
the solar limb with Hinode SOT (jTsnneta et al.l[2QM ) by 
iLites etld ] (l2nT?l revealed a thin, sub-arcsecond layer 
of photospheric emission in the 6301.5A and 6302.5A 
lines of neutral iron. These lines, normally in absorption 
within the solar disk, undergo a transition into pure emis¬ 
sion above the limb. The transition between absorption 
and emission is characterised by the appearance of emis¬ 
sive wings in the line profiles, their significant Doppler 
shifts, and linear (radially directed) light polarisation. 

Realistic numerical simulations of solar surface 
magneto-convection have long been successful in re¬ 
producing, explaining and e ven pr e dictin g various so¬ 
lar photospheric (see e.g. and refer¬ 
ences thnremJj_anTrecentlyj__chTom^ phenom¬ 

ena ([Leenaartsi^r^nilOllI; (Pereh^eT^UllOllI). Centre- 
to-limb variation of solar radiatio n has been exten¬ 
sively studied using simulations (iKeller et al.l I2QQ4 
iKoesterke et al.l I2QQ8I: iKitiashvili et al.l l2Q14l ). Ulti¬ 
mately, the success of numerical modelling in under¬ 


standing the physics of the solar photosphere can be 
demonstrated by an excellent agreement between the var¬ 
ious numerical codes used to simulate the solar surface 


In this paper, we use numerical simulati on of the so¬ 
lar p hotosphere with the MURaM code (jVogler et al.l 
l2QQ5f ) and one-di mensional line profile synth esis with the 
NICOLE code (|Socas-Navarro et al.l 12014 ) to perform 
the modelling of the spectro-polarimetric properties of 
the radiation at the solar limb with the 6301.5A and 
6302.5A Fel lines in the LTE approximation. We show 
the average properties of simulated radiation, such as 
centre-to-limb variation of the continuum intensity and 
the intensity in the line core. We also show the connec¬ 
tion between the line core intensity and the linear po¬ 
larisation above the limb. At the limb, we demonstrate 
the presence of the transition layer between pure absorp¬ 
tion and pure emissi on, spectropolarim etrically similar to 
the observations by iLites et ahl (j2QlQf ) . In the transition 
layer, the profiles with emission wings are found. Using 
response functions to the temperature for the 6301.5A 
line, we analyse the line formation process and provide a 
simple model for emission-absorption line profile forma¬ 
tion at the solar limb. In our simulations, the emission 
wings in the profiles are produced by positive tempera¬ 
ture gradients and torsional line-of-sight motions within 
the intergranular magnetic flux tube s. In contrast to the 
proposed role of non-LTE scattering (jLites et al.ll2QlQf ) in 
formation of photospheric line profiles above the limb, we 
are able to reproduce the wing emission in the profiles 
in pure one-dimensional LTE approximation. We sug¬ 
gest that both non-LTE scattering and Doppler-shifted 
emission in the flux tubes are responsible for the complex 
Stokes-/ profile shapes at the solar limb. 

The paper is organised as follows. In Section 2 we 
describe the simulation setup we used to model the solar 
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limb radiation and polarization information. Section 3 is 
divided into subsections 3.1 and 3.2, where we describe 
the general properties of the simulated limb, and analyse 
the physical origins of emission in the line profile wings at 
the limb, respectively. Section 4 concludes our findings. 

2. SIMULATION SETUP 

We use a single extended photospheric MHD model 
snapshot from a unipolar magnetic field simulation 
of a plage region produced by the MURaM code 
(jVogler et al.l l2QQ§ ). The code solves radiative MHD 
equations on a Cartesian grid. The magneto-convection 
simulation is set through injecting a uniform magnetic 
field of 200 G into a previously calculated non-magnetic 
convection model. The resolution is 25 km in the hori¬ 
zontal directions and 10 km in the vertical direction. The 
simulation box has a size of 12 x 12 Mm^ and 3.2 Mm in 
the horizontal and vertical directions, respectively. The 
average continuum formation layer is located at about 
2 Mm above the bottom boundary. The bottom bound¬ 
ary of the domain is open for in- and outflows, and the 
top boundary is open for outflows, while the inflows are 
kept at constant temperature 6000 K in order to mimick 
the chromospheric layers of the Sun. The side bound- 
aries are pe riodic. We use FreeEOS equation of state 
(jlrwinll201^ to account for changes in the gas law due to 
partial ionization. No Hall term or ambipolar diffusion 
is included in the simulation. 

A solar limb region is constructed by adding together 
16 copies of the chosen snapshot and locally changing 
the geometry of the snapshots to cylindrical geometry 
with radius Rq = 695 Mm in order to mimic the curva¬ 
ture of the solar surface in one dimension. The final size 
of the domain is about 192 x 10 Mm in the horizontal 
and vertical directions, respectively, which corresponds 
to ±8° inclination between the surface and the line-of- 
sight (LOS), or // = 0 — 0.14. The LOS components of 
the magnetic field and the plasma velocity are recalcu¬ 
lated accordingly. After the domain is projected back to 
the Cartesian grid, it is resampled to a somewhat lower 
(25 km) vertical resolution to reduce the computational 
cost, and smoothed along the lines-of-sight with a boxcar 
function with the width of 100 km to avoid numerical in¬ 
stabilities due to sharp gradients between granular and 
intergranular solar magneto-convection features during 
radiative transport computations. With the optical ray 
length of about 5 Mm the effect of the smoothing is neg¬ 
ligible. Using the same model snapshot to cover the solar 
limb allows us to study the same features under the dif¬ 
ferent inclination angles as well as at the solar limb. A 
vertical cut through half of the constructed solar limb 
domain is shown in Fig. [TJ As is evident from the fig¬ 
ure, in the constructed domain both the regions of the 
solar disk (with optically thick, hot and dense plasma) 
and the solar limb (optically thin) are present. The LOS 
components of the magnetic field and the velocity are 
predominantly horizontal in this simulation. Therefore, 
granular upflows and magnetised downflows do not in¬ 
fluence the radiation. However, the photospheric part of 
the domain is covered with horizontally expanding mag¬ 
netic fields and with the radially elongated line-of-sight 
velocity structures of alternating sign corresponding to 
Alfvenic horizontal torsional motions. These motions 
have been recently demonstrated to exist in the inter¬ 


granular ma gnetic flux con c entra ti ons in simulations and 
analys ed bv iShelvag et al ] (j2Q13f ): IShelvag fc Przvbvlskil 

(Iml . 

In order to calculate Stokes profiles for 6301.51A 
and 6302. 49 A Fel absorption lines, the horizontal one¬ 
dimensional columns from the solar limb model were used 
as an input for spectro-polarimetric line profile synthesis 
code NICOLE ([Socas-Navarro et al.lDOOOtl^cas-Navan^ 
I 2 OIII : iSocas-Navarro et al.ll2014i r~ No multi-dimensional 
scattering or non-LTE effects were included in the com¬ 
putation. The excitation potentials V = 3.65 eV and 
log gf = —0.59, and V = T686 eV and \og gf = —1.13 
were used for 6301.5lA and 6302.49A Eel lines syn¬ 
thesis, respectively. The full Stokes vectors for the 
lines were calculated within 6301.2A — 6301.8A, and 
6302. 2 A — 6302.8A wavelength ranges with 3 mA res¬ 
olution. 

3. RESULTS 
3.1. General properties 

The average radiation intensity structure across the 
simulated limb is shown in Fig. [2l The solid curve rep¬ 
resents the normalised continuum intensity variation. A 
region with sharp intensity decrease from 0.5 to nearly 
zero at 696.6 Mm has thickness of about 200 km. At 
the rest wavelength of the Eel line core, A = 6302.49A, 
the average intensity (dashed curve) is about 0.5 of the 
continuum intensity at the solar disk. However, due to 
the increased opacity in the line core, the simulated solar 
limb is higher in the solar atmosphere by about 300 km 
and spreads over a somewhat larger height range of about 
300 km. Similar behaviour is demonstrated by the min¬ 
imum intensity in the line profile (shown as dash-dotted 
curve). It should be also noted that the averaging has 
been done over a single computational snapshot, which 
does not account for the oscillations and time depen¬ 
dence present in the simulations. Therefore signatures of 
intensity variability due to granulation are also present. 

Average continuum intensites across the full solar disk 
were also computed. Fig.[3]shows the average 6300A nor¬ 
malised intensities calculated for the same model snap¬ 
shot at different observational angles at the solar disk 
// = 1, v^/2, 0.5, 0.2, 0.1, 0.05 using the domain incli¬ 
nation method (see e.g. IShelvag fc Przvbvlskil20l4 black 
solid line), the average intensities from the limb simula¬ 
tion (/i = 0.14 — 0; green solid line), the granulation 
contrast (black dashed line) and the granulation contrast 
from the limb simulation (green dashed line). As the plot 
shows, the limb darkening dependences on g fit exactly in 
the overlapping region between // = 0.05 — 0.1. The rms 
contrast dependence on g for the limb simulation shows 
large scatter due to significantly smaller statistics. How¬ 
ever, the scatter is found to be close around the average 
contrast values for the solar disk computation. 

A continuum intensity image (Fig. SI top-left panel) 
shows low-contrast solar limb granulation and a corru¬ 
gated solar limb surface. The limb height variability due 
to granulation is small, of the order of 100 — 200 km, as 
the (deeper) intergranular lanes are completely hidden by 
the yannlation. Brigh t granular walls (faculae) are visi¬ 
ble (|Keller et al.l 1200# . There are no noticeable contin¬ 
uum intensity features above the solar limb. The 6301.5A 
line core intensity image (top-right panel in Fig. S]) and. 
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Fig. 1.— Vertical cut through a half of the simulation domain representing the solar limb. The temperature, LOS velocity and LOS 
magnetic field strength are shown in the upper, middle and bottom panels, respectively. The vertical axis in the plots corresponds to the 
distance from the solar centre. The dashed line in the panels indicates approximate level of the solar surface. 
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Fig. 2.— Limb darkening across the solar limb. 6300A continuum 
intensity is shown with the solid line, the dashed line corresponds 
to the intensity at the wavelength of the line core 6302.49A, and 
the dash-dotted line corresponds to the minimum intensity of the 
Fel 6302. 49 A line profile. The horizontal axis is the distance from 
the solar centre. 

more pronounced, the 6301.5 — 0.15A line wing intensity 
image (middle-left panel), however, show a large number 
of small-scale vertically extending features with some¬ 
what enhanced line core and strongly enhanced line wing 
intensities above the average solar limb radius at this 
wavelength, R > 697 Mm. The 6302.5A line core in¬ 
tensity (middle-right panel) shows similar, although less 
prominent features. Reversed granulation is also visible 
(|Cheung et al.|[2QQ^ . 

As can be seen from the bottom panel of Fig. [U the 
LOS magnetic field mainly corresponds to the horizontal 
magnetic field of the intergranular magnetic field con¬ 
centrations, expanding in the solar atmosphere. There¬ 
fore, along the line of sight, both polarities of the mag¬ 
netic field, as well as a strong component of the magnetic 
field perpendicular to the line of sight, are present. To¬ 
gether with torsional oscillatory flows in the magnetic 


Fig. 3.— 6300A continuum limb darkening across the solar disk. 
Solid curve - simulated average normalized intensity, dashed curve 
- rms contrast of granulation. The solar limb model (/r = 0.15 — 0) 
is shown in green. Triangles show the inclination angles used to cal¬ 
culate the radiation parameters, = 1, \/3/2, 0.5, 0.2, 0.1, 0.05. 

flux tubes, this leads to generation of strong and asym¬ 
metric, multi-lobed profiles in all polarization compo¬ 
nents, as it is demonstrated in the next subsection. 

The line parameters, namely the linear 
(f + U‘^dX) and circular (/|R|(iA) polarisa¬ 
tion, are shown in the bottom-left and bottom-right 
panels of Fig. IH respectively. The images show the 
presence of strong Stokes signals both at the solar limb 
and in the simulated solar disk region. As expected 
due to the weakness of the LOS magnetic field and 
strong magnetic field component perpendicular to the 
line-of-sight above the limb, the signal in Stokes-R 
becomes weaker, and the signal in ^/Q‘^ becomes 

stronger above the solar limb region. 

Scatter plots of the intensity in the line core versus 
total Stokes-R and areas are shown in the 

left and right panels of Fig. [3 respectively. Both figures 
clearly show two separate populations of points: above 





















4 


S. Shelyag 



Fig. 4.— Continuum (top left), 6301.SA line core (top right), 6301.5 — 0.15A line wing (middle left), 6302.5A (middle right) line core 
normalised intensity images from the simulations, and f \/U‘^d\ (bottom left) and / |y|dA (bottom right) simulated images. Small- 
scale local brightenings are visible in 6301.5A line wing, and, to some extent, in the line cores. Integrated Stokes parameters show presence 
of both circular and linear polarisation. The linear polarisation is significantly more pronounced above the simulated limb compared to the 
circularly polarized light in the same region. 



Fig. 5. — Scatter plots of intensity in the line core vs the total Stokes-V area (left panel) and of intensity in the line core vs the total 
Stokes-+ U‘^ area (right panel) for Fel 6302.5A line. The green dots correspond to the r > 696.8 Mm in Fig. |4] 


the solar limb (marked by green dots) and at the so¬ 
lar disk (black dots). While the solar disk population 
does not show a dependence of the Fel 6302.5A line core 
intensity on the magnetic field, the solar limb popula¬ 
tion demonstrates such dependence: the larger the lin¬ 
ear and/or circular polarisation, the higher the line core 
intensity is. 

The horizontally averaged Stokes-/, Q//, U/I^ and 
V/I profiles for 6301.5A and 6302.5A lines across the 
solar limb are shown in Fig. [6] for the radii R = 


695.875, 696.6, 696.675, 696.775, 697.2 Mm. The pure 
absorption profiles are obtained for the position within 
the simulated solar disk at i? = 695.875 Mm. Stokes- 
Q and V show the presence of strong linear (verti¬ 
cal) and line-of-slight polarisation signals. At the limb, 
R = 696.6 Mm, Stokes-/ profiles show an absorption 
core, while the profile wings show emission. 

Stokes-y profiles have lower amplitude compared to 
the disk due to the inclination of magnetic field with 
respect to the LOS, and the multi-lobed structure of 
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Fig. 6 .— Horizontally averaged Stokes-/, Q//, t///, Vj! profiles 
for the radii R = 695.875, 696.6, 696.675, 696.775, 697.2 Mm 
(from the top to the bottom) in the model. The dotted lines show 
zero levels for Q//, t/// and Vjl profiles, which are normalised so 
one tick interval in the corrsponding plots is equal to 0.03. 

the profiles becomes more pronounced due to presence 
of both positive and negative polarities of the magnetic 
field and velocity variations along the LOS and within 
the horizontal averaging element. 

Above the limb, Stokes-/ become gradually more emis¬ 
sive in the optically thin part of the simul ated atmo- 
sphere. In agreement with the observations (|Lites et al.l 
the 6302.5A line profile becomes entirely emissive 
closer to the limb than the 6301.5A profile. At about 
300 km above the simulated limb both lines are in the 
emission regime. In the topmost layers of the simula¬ 
tion domain, both lines become fully emissive with no 
core absorption. The line intensity ratio of 3.4, which is 
close to the ratio of the values of the oscillator strengths 
used for the lines ^ 3.47 (see Section 1) indicates opti¬ 
cally thin emission regime of the lines. Stokes-Q profiles 
demonstrate strong linear polarisation above the limb, 
while Stokes-y becomes somewhat smaller, complex and 
multi-lobed, indicating that the LOS velocity and mag¬ 
netic field experience significant changes in the region of 
the line formation. Stokes-1/ does not show any signifi¬ 
cant horizontal polarisation due to the geometry of the 
magnetic field in th e simulation. 

iLites et ahl (j2QlQf ) suggested that the Fel 6301.5A and 
6302.5A line profile reversals to the emission regime in 
the wings in proximity to the solar limb are caused by 
non-LTE excitation of iron atoms by scattered light com¬ 
ing from the photosphere. Our simulation, however, 
shows the presence of such horizontally average profiles, 
computed using the pure one-dimensional LTE approx¬ 
imation. To determine the physical mechanism of the 
wing emission, we first performed the computation of the 
Ee I line profiles in the horizontally averaged model with 
zero LOS velocity. This simple numerical experiment 
showed the absence of emission wings in the computed 
line profiles and LTE-type transition of the lines from 
the absorption to emission regime across the simulated 
solar limb. Therefore, in the simulation we present here, 
the wing emission in the horizontally-averaged profiles is 
created by the superposition of the local, possibly asym¬ 
metric emission-absorption profiles produced by strong 
LOS (tangential to the solar surface) velocity and tem¬ 
perature gradients at the corrugated surface of the solar 
limb. In the next subsection, we analyse the details of 


formation of these profiles in our simulations. 

3.2. Absorption-emission profiles at the solar limb 

Eigure [7| shows simulated /, Q//, U/I and V/I slit 
images for the horizontally-averaged simulated profiles 
(top row), for a slit positioned vertically across the so¬ 
lar limb at the position marked by the vertical dotted 
line in Eig. |4] (middle row), and for the slit positioned 
horizontally at the solar limb, as marked in the top-left 
(Stokes-/) panel of the figure (bottom row). 

The horizontally-averaged Stokes-/ profiles show ab¬ 
sorption within the simulated solar disk, the emission 
profiles far above the solar limb, and the transition 
from pure absorption to pure emission, which is charac¬ 
terised by the presence of the profiles with the absorption 
core and the raised, emissive wings. Stokes-Q// profiles 
within the disk are positive predominantly because the 
magneto-convection model used in this simulation has a 
unipolar vertical magnetic field. Above the limb, Stokes- 
Q/I change their signs. Horizontally-averaged Stokes- 
U/I profiles generally have multi-lobed shapes and lower 
amplitudes, as has already been shown, due to the ab¬ 
sence of a preferred direction for the magnetic fields par¬ 
allel to the limb in the simulation. Stokes-E// profiles, 
while also being complex and multi-lobed, show some 
preference towards positive blue lobes at the disk, which 
is consistent with the geometry of magnetic field in open¬ 
ing magnetic flux tubes of positive polarity. 

Eor the vertically-aligned slit, the Stokes-/ slit image 
shows strong continuum variability on the disk and emis¬ 
sion cores above the limb for both the lines. However, 
the most pronounced feature in the image is the presence 
of bright features in the blue wings^ of both 6301.5 and 
6302.5A Ee l lines with the intensities higher and spatially 
unrelated to the continuum intensity enhancements. At 
the same locations as the bright features in the Stokes-/ 
wings, Stokes-Q// and U/I show strong polarisation sig¬ 
nals and rapidly change sign within the same intensity 
feature. Notably, Stokes-Q// and U/I are, in general, 
of similarly large amplitudes (±0.4), indicating the pres¬ 
ence of strong perpendicular to the LOS magnetic fields. 
Above the limb, Stokes-Q// show a multi-lobed struc¬ 
ture, which does not vary with height significantly, while 
Stokes-1/// vanishes. Stokes-E// signals within the disk 
show multiple lobes of both polarities too, indicating the 
presence of th e gradients of LOS v elocity and magnetic 
fields (see e.g. iShelvag et al.l [200^ . It should be noted 
that no averaging or image degradation has been applied 
for the shown vertically-aligned slit, therefore the only 
source of strongly asymmetric multi-lobed Stokes-E pro¬ 
files is the variation of the physical parameters along the 
LOS. 

Above the limb, Stokes-E// vanishes in agreement with 
the bottom-right panel of Eig. |4l 

The Stokes-/, Q//, U/I and V/I images for the 
horizontally-positioned slit at the solar limb are shown 
in the bottom row of Eig. 0 The Stokes-/ image shows 
the presence of profiles with emission wings almost ev¬ 
erywhere along the slit for both lines. It should be noted 

^ Note that the presence of only the blue wing emission here is 
purely a selection effect introduced by choosing a position with a 
bright feature in the blue wing of the profile. Similarly frequent 
features in the red wing are also found in the simulations, as well 
as the features with two emissive lobes of lower amplitude. 
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Fig. 7.— Simulated /, Q//, t///, Vjl slit images for the horizontally averaged profiles (top row), for a position across the simulated solar 
limb marked by the dotted line in Fig. 0] (middle row), and for a slit positioned horizontally at the limb (bottom row, marked in the top 
row Stokes-/ image with the dashed line). The position of the vertical slit shown in the middle row is marked by the horizontal dashed line 
in the bottom row. 


that most of the profiles are strongly asymmetric, with 
emission either in the blue or in the red lobe. The pres¬ 
ence of the emission lobes shows some correlation with 
the continuum intensity variation. The emission lobes 
are strongly shifted with the corresponding LOS veloci¬ 
ties of 4 — 8 km s“^. There is a weak visual correlation 
between Stokes-Q//, Ujl amplitudes and the emission 
lobes, while Stokes-Q// has expectedly larger amplitude 
than Stokes-t///. Nevertheless, the Stokes-Q// and Ujl 
signal patterns coincide with the emission lobes of the 
profiles, while there are no strong signals in the line ab¬ 
sorption cores. A similar dependence is found for the 
Stokes-f/// signal. 

The process of formation of the profiles which show 


both emission and absorption is analysed below. An 
arbitrary position within the simulated solar disk, but 
close to the limb was selected, where the 6301A line 
profile shows both absorption in the line core and emis¬ 
sion in the blue wing. Response functions (RFs) of all 
Stokes parameters to a small temperature perturbation 
(AT = 100 K) in the corresponding one-dimensional at¬ 
mosphere were computed. Fig. [8] shows the calculated 
RFs (grayscale images) together with the correspond¬ 
ing line profiles (plots in the top row), and the photo- 
spheric parameters (vertical plots). A clear correspon¬ 
dence is seen between the emission part of the profile 
at A A = —0.1 A, and a sharp temperature rise at about 
1.5 Mm distance along the LOS. At the same location. 
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Fig. 8 .— Simulated /, Q//, U/I^ V/I profiles (top row), their response functions to the temperature perturbation (grayscale images) and 
the structure of the corresponding ID atmosphere along the LOS. Dashed vertical lines mark zero values for vios^ Bios^ and B±. 


a strong plasma flow of about 4 — 5 km s“^ towards the 
observer is found, which shifts the emission wavelength 
towards the blue side of the spectrum. The LOS com¬ 
ponent of the magnetic field Bios changes the direction 
and amplitude from —200 G to 250 G, leading to the 
formation of a strong and asymm etric Stokes-V7/ sig - 
nal at the same wavelength (see e.g. lShelvag et al.ll2007[ ). 
The vertical, perpendicular to the LOS component of the 
magnetic field is almost always positive with the min¬ 
imum magnetic field of a few Gauss. In the emission 
region it increases up to about 600 G, indicating that 
the line-of-sight crosses an intergranular magnetic field 
concentration there. The Stokes-I/// RF has a maxi¬ 
mum in this region, and the Stokes-I/// profile shows 
two lobes of large amplitude in the blue. The horizon¬ 
tal, perpendicular to the LOS magnetic field component 
B± has expectedly smaller amplitude and also shows an 
increase in the emission region up to 250 G. 

A summary of the process of photospheric emission- 
absorption Fe I line profile formation is provided in Fig.O 
The line-of-sight crosses three distinct regions. The first 
region (A) is above the grannie, where both the LOS 
components of the veiocity and the magnetic fieid are 
smaii (from 0 to about 1.0 Mm aiong the LOS in Fig. [51 
In this region, the continuum radiation is formed. Then, 
the line-of-sight crosses the intergranular magnetic field 
concentration region (B). Strong torsional oscillatory mo¬ 
tions with spe eds of a few km s~^ within the magnetic 
conce ntration (jShelvag et al.ll20T3l: iShelvag fc Przvbvlskil 
l2014f ) Doppler-shift the line formation wavelength in ac¬ 
cordance to the direction of motion in this region (from 



Fig. 9. — Emission-absorption line profile formation in the non- 
uniform solar photosphere at the limb. The line-of-sight close to the 
solar limb crosses three different regions: (A) a non-magnetic gran¬ 
ular region, where the continuum radiation is formed, (B) the inter¬ 
granular magnetic field concentration, where the emission part of 
the profile is formed, and which is characterised by radiative heat¬ 
ing, positive temperature gradient (T 2 > Ti), and torsional plasma 
motions (id), and (C) a non-magnetic region over the granule with 
a negative temperature gradient (T 3 < T 2 ), where the absorption 
part of the profile is formed. 

1.0 to 2.0 Mm in Fig. [8]). At the same time, radia- 
tive heat exchange with the hot granular wall (see e.g. 
iSprnitI Il976l : iVogler et ah! l2004f ) and resistive heating 
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(|Mq11 et al.ll2Q12l: iShelvag fc Przvbvlskill2Q14l) cause the 
localised temperature increase within the magnetic flux 
tube and emission. Finally, the line-of-sight leaves the 
magnetic field concentration and enters the rarified re¬ 
gion above the next granule (C), where the LOS veloci¬ 
ties are generally smaller and temperature decreases (up 
to the chromospheric layers, to which the photospheric 
lines under consideration are not sensitive), leading to 
the formation of an absorption line core. 

4. DISCUSSION 

In this paper, we presented a novel spectro- 
polarimetric simulation of photospheric radiation at the 
solar limb. We used a single snapshot from magneto¬ 
convection simulations of a plage region to construct a 
cylindrical solar limb-like region and performed detailed 
radiative diagnostics computation on this domain with 
the 6301.5A and 6302.5A pair of photospheric lines of 
neutral iron in the LTE approximation. We computed 
the average radiation parameters at the transition be¬ 
tween the solar disk and limb, such as limb darkening, 
rms granulation contrast and average Stokes profiles. 

The average Stokes profiles showed the emission above 
the limb, the presence of emissive line wings, absorption 
core and significant pola risation signal at the limb, sim¬ 
ilar to the observations (|Lites et al.ll2QlQ[ ). where it was 
suggested that the emissive wings are the result of ra¬ 
diation scattering and non-LTE excitation in the solar 
atmosphere. Our results, however, suggest that the ra¬ 
diation scattering at the limb is not the only mechanism 
causing the emission in the line wings. It is demonstrated 
that, in the LTE approximation, the emission lobes in 
Stokes-/ profiles at the limb can be formed due to the 
presence of the torsional bidirectional flows along the 
line-of-sight and the local heating in the intergranular 
magnetic field concentrations. When the line-of-sight, 
along which the radiation is formed, crosses the inter¬ 
granular magnetic concentration region, it experiences a 
temperature rise, leading to emission, and a Doppler shift 
to the blue (or red) side of the spectrum. After that, the 
line-of-sight crosses an absorption region with a smaller 
LOS velocity and a temperature decrease, which leads 
to the formation of the absorption core. As there is no 
preferred direction for the torsional flows in magnetic 
field concentrations, both blue- and red-wing emission 
is observed in the simulations. The profiles with two 
emissive lobes are also found in the regions where the 
line-of-sight crosses t he layers of opposite r otation in the 
magnetic flux tubes (|Shelvag et al.ll2QllbD . When aver¬ 
aged spatially, these profiles lead to the appearance of the 
Stokes-/ profiles with emissive lobes and absorption core. 
The polarisation signals at the limb show the presence of 
vertically (radially) oriented magnetic fields, which is ex¬ 
pected since a plage model with 200 G vertical magnetic 
field was used for the computation. 

Despite the detailed characteristics of torsional mo¬ 
tions in magnetic field concentrations still need to be 


better understood (see e.g. IWedemever fc SteineJl2014D . 
there is little doubt of their presence in the mag¬ 
netic field concentrations. Eurthermore, non-magnetic 
or weakly magnetised photospheric mod els also show 
torsional motions in int e rgranular lanes (l^elvag et ^ 
I2011al IMoIT et al ] I2Q11I: iKitiashvili et al.l I2Q12L ]2Q13f ). 
While in this case of zero or weak photospheric magnetic 
field the magnetic tension is small and does not prevent 
tornado-like appearance of the vortex, the purely hydro- 
dynamic vortex flows in the intergranular lanes would 
lead to the presence of torsional motions consequently 
leading to the similar mechanism of emission-absorption 
line profile formation at the solar limb. 

A natural point of concern is the use of LTE approxi¬ 
mation in the presented study. While there is no doubt 
the non-LTE effects are important especially at the limb, 
the primary aim of the study was to demonstrate the LTE 
effects of the limb radiation formation. Eurthermore, the 
neutral iron lines used for spectropolarimetric diagnos¬ 
tics in this paper are formed low in the photosphere ir¬ 
respectively of the position at the solar disk, where the 
LTE assumptio n is still reasonable, at least for the inter¬ 
granular lanes (jShchnkina fc Trniillo BnenolDOOil ). The 
basic process of emission-absorption line profiles forma¬ 
tion at the limb, which is described in the paper, relies 
on the presence of torsional motions and heating within 
the intergranular magnetic flux tubes. 

The structures described in the paper are rather small- 
scale, with linear dimensions of about 100 km. Taking 
into account low light and high noise levels, inevitably 
associated with limb observations, it would be difficult 
to provide detailed and direct observational confirmation 
of the effect of torsional flows in magnetic field concen¬ 
trations on the photospheric line formation at the solar 
limb. Nevertheless, future, large-scale instruments, such 
as DKIST, will be able to observe these small features. 
Einally, further study is needed and planned to analyse 
their lifetimes and links to chromos pheric fine structu re, 
bidirectional flows and oscillations (jSekse et al.ll2Q13[ ). 
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